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ABSTRACT 

Aims. We aim at deriving the physical conditions in the neutral gas associated with damped Lyman-a- systems using observation and 
analysis of IF and C i absorptions. 

Methods. We obtained a high-resolution VLT-UVES spectrum of the quasar Q 2348-01 1 over a wavelength range that covers most of 
the prominent metal and molecular absorption lines from the log ^(H i) = 20.50±0. 10 damped Lyman-o- system at z a bs = 2.4263. We 
detected H 2 in this system and measured column densities of H 2 , C i, C i*, Ci", Sin, Pn, S n, Fen, and Nin. From the column density 
ratios and, in particular, the relative populations of IF rotational and C i fine-structure levels, we derived the physical conditions in 
the gas (relative abundances, dust-depletion, particle density, kinetic temperature, and ionising flux) and discuss physical conditions 
in the neutral phase. 

Results. Molecular hydrogen was detected in seven components in the first four rotational levels (J = 0-3) of the vibrational ground 
state. Absorption lines of IF J = 4 (resp. J = 5) rotational levels are detected in six (resp. two) of these components. This leads to 
a total molecular fraction of log/ -1.69*qjI. Fourteen components are needed to reproduce the metal-line profiles. The overall 
metallicity is found to be -0.80, -0.62, — 1 . 17±0. 10 for, respectively, [Si/H], [S/H] and [Fe/H]. We confirm the earlier findings that 
there is a correlation between log iV(Fe n)/N(S n) and log N(Si n)/W(S n) from different components indicative of a dust-depletion 
pattern. Surprisingly, however, the depletion of metals onto dust in the IF components is not large in this system: [Fe/S] = -0.8 to 
-0.1. 

The gas in H 2 -bearing components is found to be cold but still hotter than similar gas in our Galaxy (7* > 130 K, instead of typically 
80 K) and dense (n ~ 100 - 200 cm -3 ). There is an anti-correlation (R = -0.97) between the logarithm of the photo-absorption rate, 
logy8 , and logA , (H 2 )/A'(Ci) derived for each H 2 component. We show that this is mostly due to shielding effects and imply that the 
photo-absorption rate /3q is a good indicator of the physical conditions in the gas. We find that the gas is immersed in an intense UV 
field, about one order of magnitude higher than in the solar vicinity. These results suggest that the gas in H 2 -bearing DFAs is clumpy, 
and star-formation occurs in the associated object. 
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1. Introduction 

High-redshift damped Lyman-a systems (DLAs) detected in ab- 
sorption in QSO spectra are characterised by their large neutral 
hydrogen column densities, N(Hi) > 2 x 1 2Q cm" 2 , similar t o 
what is measured through local spiral disks (IWolfe et al.l fl986). 
The corresponding absorbing clouds are believed t o be the reser- 
voir of neutral hydrogen in the Universe (e.g., iPeroux et al.l 
l2Q03t iProchaska et al.ll200 5). Though observational studies of 
DLAs have been pursued over more than two decades (see 
IWolfe et aill2005l for a recent review), important questions are 
still unanswered, such as (i) the amount of in-situ star-formation 
activity in DLAs, (ii) the connection between observed abun- 



* Based on observations carried out at the European Southern 
Observatory (ESO) under prog. ID No. 072.A-0346 with the UVES 
spectrograph installed at the Very Large Telescope (VLT) Unit 2, 
Kueyen, on Cerro Paranal, Chile. 



dance ratios and the dust content, and (iii) how severe is the bias 
due to dust obscuration in current DLA samples. 

One way to tackle these questions is to derive the physical 
conditions that prevail in the absorbing gas by studying the exci- 
tation of molecular and atomic species together. Since UV-lines 
are redshifted into the optical, a single VLT-UVES spectrum al- 
lows one to simultaneously observe a large number of absorption 
lines of important species at high spectral resolution and with a 
good signal-to-noise ratio. In the case of the Galactic interstellar 
medium (ISM), only a few lines of sight have been observed in 
such detail. 

Although molecular hydrogen is the most abundant molecule 
in the Universe, it is very difficult to detect directly. The 
molecule has a rich absorption spectrum in the UV, and space- 
crafts have been used to detect the absorption lines from clouds 
located in our Galaxy in front of bright background stars. At 
Hi column densities as high as those measured in DLAs or 
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sub-DLAs (N(Hi) > 19.5 cm -2 ), H2 molecules are conspic- 
uous in our Galaxy: sigh tlines with logjV(Hi) > 21 usually 
have log N(U 2 ) > 19 (see lSavage et aljfl977T: iJenkins & Shaval 
11979b . More recently, observations with the Far Ultraviolet 
Spectroscopic Explorer (FUSE) have derived the physical con- 
ditions in the diffus e molecular gas in the Galaxy but also the 
Magellanic clouds (Tumlinson et al. 2002) or high Galactic lat- 
itude lines of sight dRichter et al.l 12003b IWakkerll2006l) . In the 
Magellanic clouds, molecular fractions and kinetic temperatures 
are reported to be, respectively, lower and higher than in the local 
ISM. This is probably a consequence of higher UV ambient flux 
and lower dust content. In intermediate velocity clouds located 
in the Galactic halo, H2 is found to arise in cold (T $ 140 K), 
dense («h ~ 30 c m" 3 ), and clumpy (D ~ 0.1 pc) clouds 
(Richteretal]|2003), embedded in a relatively low UV flux. 

At high redshift, UV-lines are redshifted in the opti- 
cal so ground-based t el escopes can be used. Quasars (e.g. 
Levshakov et al. Il992t iGe & Bechtokl fl999; Ledoux et al 



2001 " but also recently y-ray burst afterglows dFvnbo et al 
20061) are used as background sources. The amount of H 



molecules s een in DLAs is much less than what is observed in 
our Galaxy. ILedoux et al.l d2003l) find that about 15% of DLAs 
show detectable H2 with molecular fraction in the range -4 < 
log / < — 1. This fractio n increases with met allicity and reaches 
50% for [X/H] > -0.7 dPetitiean et alj[2006h . The fact that H 2 
is less conspicuous at high redshift can be explained as a conse- 
quence of lower metal content and therefore lower dust content, 
together with higher ambient UV flux. It could be also that the 
covering factor of the diffuse molecular gas is low. 

Detaile d analysis of H2- bearing DLAs has been pursued in 
some cases. ICui et al.l d2005l) recently derived a surprisingly low 
hydrogen density, «h ~ 0.2 cm 4 , using R ~ 30,000 STIS data 
from the HST, and an electron temperature of T e w 140 K 
in the H 2 system at z abs = 1.776 toward Q 1331 + 170. Other 
authors, however, usually derive higher densities. By studying 
the spectrum of the z a b s ^ 2 .811 DLA toward PKS 0528-250, 
ISrianand & Petitieanl ([1998) derived a kinetic temperature of 
T ~ 200 K and a density of n ~ 1000 cm 4 . They suggest that 
the ratio N(H2)/N(C 1) could be a good ind icator of the physical 
conditions in the gas. Petitiea n"et al.l(l2002l) derived gas pressure, 
particle density, excitation temperatures, and ambient UV fluxes 
for four H2 components at z a bs - 1.973 toward Q 0013-004. 
They show that, whenever H2 is detected, the particle density is 
high («h ~ 30 - 200 cm 4 ) and the kinetic temperature is low 
{T ~ 100 K). In addition, the ambient UV radiatio n field is found 
to be highly inhomogeneous. ILedoux et al.l (120021) find high par- 
ticle densities (n ~ 30 - 400 cm 4 ) and moderate UV radiation 
field in a DLA at z a bs = 1 .962. They conclude that physical con- 
ditions play an important rol e in go verning the molecular frac- 
tion in DLAs. Rei mers et al.l d2003l) also found a high particle 
density (n ~ 100 cm 4 ) in a DLA at z a b s = 1.150. From the high 
derived UV radiation field, they suggest there i s ongoing star- 
forma tion activity close to the absorbing cloud. Srianan d et al.l 
(2005) use the 33 systems from t he hig h-resolution (R ~ 45,000) 
UVES survey by ILedoux et al.1 d2003) - among which 8 show 
detectable H2 absorption lines - to compare observations and 
molecular cloud models. They conclude that the mean kinetic 
temperature in H2-bearing DLAs is higher (~ 150 K) than in the 
local ISM (~ 80 K), the UV flux is close to or higher than the 
mean UV flux in the ISM of the Galaxy, and the particle den- 
sities (~ 10-200 cm 4 ) and photo-ionisation rate are similar to 
the ones in the Galactic cold neutral medium. More recently, gas 
temperature, T ~ 90 - 180 K, and density, «h ^ 50 cm 4 , were 
derived in a DLA at z a bs - 4.224, which corresponds to the high- 



est redshift at which this kind of study has been done up to now 
dLedoux et al.ll2006bl) . 

We recently reported the detection of H2 in the z a b s = 2.4263 
DLA toward Q 2348-011 dPetitiean et al.l 120061) . We present 
here a detailed analysis of the physical conditions in the gas. This 
system has one of the highest molecular fractions ever measured 
at high redshift, and seven H2 components are detected in which 
physical conditions can be derived. This is the first time that a 
comparison of physical conditions can be made in such a large 
number of components in a single absorption system. The paper 
is organised as follows. In Sect. 2, we describe observations and 
data reduction. In Sect. 3, we give details on the data and present 
the fits to the H2, C 1 and low-ionisation metal absorption lines. 
In Sect. 4, we discuss the relative populations of H2 rotational 
levels, and C 1 fine-structure levels, deriving information on the 
physical state of the gas. Finally, we discuss our findings and 
draw our conclusions in Sect. 5. 

2. Observations 

The quasar Q 2348-011 was observed with the Ultraviole t 
and Visible Echelle Spectrograph (UVES; iDekker et all feOOO) 
mounted on the ESO Kueyen VLT-UT2 8.2 m telescope on Cerro 
Paranal in Chile in visitor mode on October 29 and 30, 2003, 
under good seeing conditions (FWHM < 0.8 arcsec) and at low 
airmass {AM < 1.2). Three 5400 s exposures were taken with 
Dichroic #1 simultaneously with the blue and red arms. We used 
a 1 .0 arcsec slit width and a pixel binning of 2 x 2, leading to 
a resolving power o f R g 47,000. The d ata were reduced using 
the UVES pipeline (Bal lester et al.ll2000l) . which is available in 
the context of the ESO MIDAS data reduction system. The main 
characteristics of the pipeline are to perform a precise inter-order 
background subtraction, especially for master flat fields, and to 
allow for an optimal extraction of the object signal rejecting cos- 
mic rays and performing sky subtraction at the same time. The 
pipeline products were checked step-by-step. 

The wavelength scale of the reduced spectra was then con- 
verted to vacuum-heliocentric values, and the portions of the 
spectrum corresponding to different arms were rebinned to a 
constant wavelength step. Individual 1-D exposures were scaled, 
weighted, and combined. No other rebin was performed during 
the analysis of the whole spectrum. Standard Voigt-profile fitting 
methods were used for the analysis to determine column densi- 
ties using the wavelengths and oscillator strengths compiled in 
Table 1 of ILedoux et alj (2003) f or metal ions and the oscilla- 
tor strengths given by iMorton & Dinersteinl (19761) f or H2.We 
adopted the Solar abund ances given by Morton (12003) based on 
the meteoritic data from lGrevesse & Sauvalld2002l) . 

3. The DLA at z abs = 2.4263 toward Q 2348-011 

There are two DLAs atz a b s = 2.42630 andz a b s = 2.61473 toward 
Q 2348-011. These DLAs were discovered as strong Lyman- 
a absorptions in a l ow spectral-resolution (~ 5 A) spectrum by 
Sarg ent et al.l d!989l) . who also observed strong metal lines as- 
sociated to the z a b s = 2.4263 DLA. From Voigt-profile fitting to 
the damped Lyman-o- line, we find that the total neutral hydro- 
gen column density of the system at z a bs = 2.4263 is log N(H i) = 
20.50+0.10 dLedoux et al.l2006al) . A simultaneous Voigt-profile 
fitting of the Lyman-o' and Lyman-/? lines was also performed 
for the system at z a bs = 2.6147. Note that the z a b s = 2.6147 
Lyman-/? line falls in the wavelength range where FL lines at 
z a b s = 2.4263 are detected (see Fig. [TJ. The total neutral hydro- 
gen column density for this system is log A^(Hi) = 21.30 + 0.08 
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dLedoux et al.l2006al) . No molecular hydrogen is detected at this 
redshift down to log 7V(H2) < 13.7 (considering rotational levels 
J = 0andJ= 1), or log/ < -7.2. 

3.1. Molecular content 

The H2 absorption lines at z a b s = 2.4263 are strong and numer- 
ous (see Fig.[TJ. Although they are all redshifted in the Lyman-o 
forest, the resolving power of our data (R st 47,000) is high 
enough to allow unambiguous detection and accurate determi- 
nation of the line parameters. 

3.1.1. Overall Fit 

The absorption profiles of the singly-ionised metal species 
have been fitted using fourteen components and, in the follow- 
ing, we use their numbering to refer to any individual veloc- 
ity component. Seven of these components, spread over about 
300kms _I , show associated H2 absorption. Multiple component 
Voigt-profile fittings have been performed simultaneously on the 
H2 absorption lines seen redward of the Lyman break produced 
by the z a b s = 2.6147 system, also including the Voigt profile of 
the Ly-y6 line at z a b s =2.6147. The continuum normalisation was 
performed using first-order broken lines between 3600 A 

and /l bs - 3850 A and a routine based on the IRAF task contin- 
uum that can fit segments of 50-100 A at once. For each ~ 50 A 
wavelength range, the continuum is fitted by approximately five 
to seven joint linear segments adjusted by an iterative procedure 
that rejects absorptions over a certain threshold related to the 
local signal-to-noise ratio. This allows us to estimate the contin- 
uum position where the Lyman-o forest is dense, and over the 
Zabs = 2.6147 Ly-/? absorption (see Fig.[T). Note that only the H2 
Lyman band is covered by our spectrum. 

The solution of a multi-component Voigt-profile fit is often 
not unique, especially when absorption lines are blended and 
sometimes saturated. In that case, it is difficult to constrain both 
the Doppler parameters and the column densities at the same 
time for all species. We therefore obtained a first guess of the 
molecular parameters by fitting selected lines (see Table. [TJ of 
the first four rotational levels, J = to J = 3, using Doppler pa- 
rameters and redshifts derived from the C 1 profiles, which are 
not saturated and located outside the Lyman-o- forest. We then 
froze the corresponding parameters and added the J = 4 and J = 5 
features. After a satisfactory fit to all H2 transitions had been 
found in these conditions, we relaxed all parameters one by one, 
until the best fit was obtained. Then, a range of column densities 
was derived by fitting the absorption features with the extreme 
values of the Doppler parameter (see Table The errors given 
in column 4 of Table [2] therefore correspond to a range of col- 
umn densities and not to the RMS errors derived from fitting 
the Voigt profiles. The quality of the overall fit can be assessed 
from Fig. [T]and in Figs.|9]to[T4] Results are presented in Table|2] 
The total H2 column density integrated over all rotational levels 
is log MFL;) = 18.52^-^, corresponding to a molecular fraction 
log / = -1.69*9 5Z. We also derive a 3 cr upper limit on the col- 
umn density of HD molecules, leading to log 7V(HD)//V(H2) < 
-3.3 for the strongest H2 components (# 1, 9 and 10). We de- 
scribe individual components below. 

3.1.2. Comments on individual components 

Components 1 & 2, z abs = 2.42449(7) and 2.42463(8): These two 
components are heavily blended and appear as a single saturated 



feature for low H2 rotational levels, but are better separated for 
J > 4, in C 1 and low-ionisation metal lines (see Fig.[2l[3land[T3l. 
Doppler parameters, b, are mainly constrained by the H2 L0- 
0R0, H 2 L1-0P2, H 2 L1-0P3 profiles, and the transitions from 
the J > 4 rotational levels, but they remain uncertain. No damp- 
ing wing is present. The first guess to fit the H2 profiles, using 
redshifts from C 1, C 1*, and C 1** transitions, was extremely help- 
ful in this case. It must be emphasised that the blends of the two 
components can be fitted using either low b values, leading to 
high column densities or higher b values, leading to lower col- 
umn densities. The best fit is obtained, however, with b = 2.9 
kms^ 1 and column densities log 7V(H 2 , J) = 17.15; 17.70; 17.00; 
17.60; 15.40 and 15.10 for J = to J = 5, respectively. These are 
the values we use in the following. Rotational levels are unam- 
biguously detected up to J = 5. Possible absorptions from J = 6 
and J = 7 levels, unfortunately, are always redshifted on top of 
strong H 1 Lyman-o- features or in regions of low SNR. Upper 
limits are given in Table [2] 

Component 3, z a b s = 2.42491(7): This component is the weak- 
est in the system and is not saturated. The Doppler parameter 
is therefore mainly constrained by the strongest transition lines 
that are not blended: H 2 L2-0 R0, H 2 Ll-0 R0, and H 2 L3-0 PI . 
Only the first four rotational levels are detected. J = 3 lines are 
always blended. However, the column-density can still be esti- 
mated using H 2 L1-0R3 (left-half of the profile), H 2 L1-0P3, 
and H 2 LO-0 R3 (see Fig.EJl. 

Components 9 & 10, z abs = 2.42659(1) and 2.42688(5): These 
components have high H2 column densities. They are separated 
well and Doppler parameters are constrained well by the width 
of the isolated and non-saturated absorption features H2 L0- 
0R0, H 2 L0-0R1 and mostly of the P-branch (AJ = -1) Lyman- 
band transitions from the J = 2 to J — 4 rotational levels. Note 
that component 9 of a particular transition is often blended with 
components 1 and 2 of another transition from the next rotational 
level (see e.g., H2L2-ORO in Fig. [9]). The first five rotational 
levels (J = to J = 4) are detected. Some consistent features are 
seen for J = 5 (see transitions H2 L4-0 P5 and H2 L2-0 R5 on 
Fig. [T4b . Upper-limits were derived as the largest value consis- 
tent with the observed spectrum for levels J = 5, J = 6, and J = 7. 
Components 12 & 13, z abs = 2.42713(5) and 2.42729(3): 
Absorptions from rotational levels J = to J = 4 are detected in 
the two components. 7V(J = 2) is constrained mainly by H2 L0- 
P2, but the maximum TV-value is also constrained by H2 L3- 
0R2 and H2LI-OP2, whereas the minimum TV- value is con- 
strained by the local depth of the features clearly seen in H2 Ll- 
0R2 and H2 L0-0R2 (see Fig. ITTb. Similar constraints are used 
for J = 3. J = 4 lines are blended or located in poor SNR re- 
gions, but we can still constrain the column-density thanks to 
H 2 L4-0R4, H 2 L4-0P4 and H 2 L3-0P4. The J = 5 rotational 
levels may be detected; but the corresponding column densities 
are very uncertain, and one should probably consider these val- 
ues only as upper limits (see Table|5J. 

A summary of the transitions used to constrain the column- 
density for each ^-component in a given rotational level is 
given in TableQ] The lines in this table are non-blended and have 
a rather good SNR, except when marked by "u" or "1" indicat- 
ing that the line is used to estimate, respectively, an "upper" or a 
"lower" limit. 

3.2. Neutral carbon 

Neutral carbon is detected in nine components. The corre- 
sponding absorption profiles are complex, as resulting from 
the blend of absorption lines from three fine-structure lev- 
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Fig. 1. A portion of the normalised spectrum of Q 2348-01 1 . The resolving power is about R 47,000. The H2 spectrum derived 
from Voigt-profile fitting is superimposed to the quasar spectrum. The strong feature at A 3707 A is the Ly-{3 absorption line from 
the DLA at z a b s = 2.6147. H2 components are numbered for the J = rotational level only. Transition names are abbreviated as 
LxYz standing for H2 Lx-OYz, and only marked for component 1. The best fit parameters are given in Table [2] 



els of the neutral carbon ground state (2s 2 2p 2 3 Po,i,2)- In 
Fig. |2] we show the fit to the C 1 absorption features around 
/lobs ~ 5346 A and A bs ~ 5676 A with contributions from the 
eleven transitions: C 1/11560, C 1/11656, Cf/LU560.6,1560.7, 
Cf/LU1656, 1657.3, 1657.9, Ci**/U1561.3,1561.4, and 
C i**/l/11657,1658. We first performed a fit with seven Ci 
components corresponding to the seven H2 detected compo- 
nents, but we left the exact positions of the C 1 components free. 
We then added Ci*, and finally Ci** transitions. To reproduce 
the overall profile well, we had to introduce two weak Ci 
components (# 6 and 7) that are not detected in H2, but are 
detected in the absorption profiles of low-ionisation metals. 



We then performed a second fit with nine components, in the 
same order: first Ci only, then adding Ci* and finally Ci**. 
Finally all transitions were fitted together. The resulting small 
velocity shifts between C 1 and H2 components (see Table [3) are 
within uncertainties. This means that H2 and C 1 absorptions are 
generally produced in the same regions. Note, however, that the 
Doppler parameters found for the C 1 components are slightly 
larger than those found for the H2 components, indicating 
that the C 1 regions are probably of slightly larger dimensions. 
Results of the fit (Fig. |2]l are presented in Table [3] Three cr 
upper-limits are derived for N(Cf*) in components 6, 7, and 9, 
leading to upper-limits on N(C f*)/N(C 1) (see Table|3]l. 
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# 


Zabs 


b 


logA'(Ci) 


logMCi*) 


log _/V(Ci") 1 


Av C i/h 2 ' 


log A^Ci'VNtC i) 


logA'tC i**)/N(C i) 






[kms '] 










n — 1 i 

[kms '] 






1 


2.42450(5) 


4.2+ 0.3 


13.55+ 0.02 


13.51± 


0.03 


13.10± 0.04 


+0.7 


-0.04±0.04 


-0.45+0.04 


2 


2.42465(4) 


4.6+ 0.4 


13.10+ 0.06 


13.06± 


0.02 


1 O C A i A 1 A 

12.54± 0.10 


+ 1.4 


-0.04±0.06 


-0.56+0.12 


3 


2.42494(4) 


2.8+ 0.4 


12.55± 0.08 


12.66± 


0.09 


12.60± 0.12 


+2.4 


0.11+0.12 


0.05±0.14 


6 


2.42605(3) 


5.9± 0.5 


12.70+ 0.07 


12.50± 


0.09 


<12.50 




-0.20±0.11 


<-0.20 


7 


2.42628(9) 


8.4± 0.5 


12.96± 0.04 


12.75± 


0.04 


<12.50 




-0.21±0.06 


<-0.46 


9 


2.42660(0) 


4.0± 0.5 


12.80± 0.02 


12.93± 


0.10 


<12.50 


+0.8 


0.13+0.10 


<-0.30 


10 


2.42688(9) 


4.0± 0.5 


13.09± 0.02 


13.05+ 


0.04 


12.55+: 0.20 


+0.3 


-0.04±0.04 


-0.54±0.20 


12 


2.42715(1) 


3.2± 0.5 


12.87+ 0.04 


12.84± 


0.04 


12.53+ 0.05 


+ 1.4 


-0.03±0.06 


-0.34±0.06 


13 


2.42727(3) 


3.0+ 0.6 


12.80+ 0.05 


12.90± 


0.03 


12.80± 0.10 


-1.7 


0.10+0.06 


0.00+0.11 



1 3 tr upper-limits are derived for A^C i**) in components 6, 7 and 9. 

2 We estimate the errors on Avci/Ht to be about 2 to 3 kms -1 . 

3 There is no velocity shift between the components of neutral carbon and those of molecular hydrogen. 



Table 1. Transitions used to constrain H2 column densities in 
the seven components. 



Level 


#1 &2 


#3 


#9& 10 


#12 


#13 


J = 


L0R0 




L0R0 


L0R0 


L0R0 




L1R0 


L1R0 


L1R0 


L1R0 






L2R0 


L2R0 


L2R0 


L2R0 


L2R0 


J = 1 


L0P1 


L0P1 










L0R1 


L0R1 


L0R1 


L0R1 


L0R1 




L1P1 


L1P1 


L1P1 








L1R1 




L1R1 


L1R1 






L2R1 




L2R1 


L2R1 






L3P1 


L3P1 














L3R1 


L3R1 




J = 2 


L0P2 




L0P2 


L0P2 


L0P2 




L0R2 


L0R2 




L0R2' 


L0R2' 




L1P2 


L1P2 


L1P2 


L1P2 


L1P2" 




L1R2 






L1R2' 


L1R2' 




L3P2 


L3P2 
















L3R2" 


L3R2" 


J = 3 




L0R3 
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L1R4" 










L2P4 
















L3P4 


L3P4 


L3P4" 




L4P4 
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L4P4 


L4P4 
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L4R4 


L4R4" 




L4R4 


L4R4 


J = 5 
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L2R5" 


L2R5" 






L3R5 


L3R5" 














L4P5" 


L4P5" 


L4P5" 



Transition names are abbreviated as LxYz standing for H2 Lx-0 Yz. 
"u" (resp. "1") means that the transition is used to give an upper (resp. 
lower) limit on N(H 2 ; J). 

3.3. Metal content 

The absorption profiles of low-ionisation species (P n, S 11, Fe 11, 
Nin and Sin) are complex and span about Av ^ 300 kms -1 in 
velocity space. Fourteen components were needed to reconstruct 
the whole profiles, i.e., the nine C 1 components plus five com- 
ponents that are not detected in C 1. A few profiles are shown in 




123 67 9101213. 



5342 5344 5346 5348 5350 5352 

Observed wavelength [A] 




12 3 6 7 9 101213 

III I I I I I 

123 679101 21 3 

III I I I I I I 

12 3 6 7 9 101213 



'. 1 2 3 6 7 9 101213 ~ 

5672 5674 5676 5678 5680 5682 5684 
Observed wavelength [A] 

Fig. 2. Observed portions of the normalised spectrum of 
Q 2348-011 (upper part in both panels) where Ci absorptions 
are seen, with the model spectrum superimposed. The modelled 
spectrum is decomposed into C 1 (solid line), C 1* (dashed line), 
and C 1** (dotted line) model spectra (lower part in both panels). 
The positions of the nine components are marked by small ver- 
tical lines. The associated labels correspond to the numbering in 
Table[3](see also Tables |2]and|4]i. 



Fig. [3] and results from the fit are given in Table [4] The overall 
metallicity is found to be -0.80, -0.62, and -1.17+0.10for, re- 
spectively, [Si/H], [S/H], and [Fe/H]. The strongest component 
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Table 2. Voigt-profile fitting results for different rotational lev- 
els of the vibrational ground-state transition lines of H2. 



# 



-Cabs 



[kms" 1 ] 



"Level log _/V(H 2 , J) 2 



b 

[km s 



1 2.42449(7) 
-158 



2 2.42463(8) 
-145 



3 2.42491(7) 
-121 



9 2.42659(1) 
+25 



10 2.42688(5) 
+51 



12 2.42713(5) 
+73 



13 2.42729(3) 
+87 



1:20 
-0.20 



15.55-17.40 
16.40-18.00 
15.90-17.70 
16.10-17.80 
15.25-15.50 
14.90-15.20 

< 14.00 
<13.80 
15 35 +01 ° 

15 60 + ™ 
15 30 +!l3!i 
14.55 + °!° 
14.30+ 
<13.80 
<13.70 

14 60 +0 10 

15 00 +0 05 

J -0,50 

14 30 +M ° 
14 55+0 ™ 

<13.80 
<13.70 
17 30 +0 2 ° 
1 7 7O +0 10 

16 50 +a9 ° 

16 70 +075 
14 55+0-30 

< 14.00 

< 13.50 
<13.70 

17 30 +0 - 20 

-P.6P. 
17 90 +il2 <> 

17 40+° 20 

17.00+g 

14.85+g 

< 14.00 

< 13.50 

<13.70 

14 95 +o.io 

15 35 +m 

14 QQ+O tO 

15 10 +0 10 

14 20 +0 10 
it.z.u _ 030 

<13.70 

14 65 +0 05 

15 35 M ° 
15 00 +0 10 
15 20 + ° 10 
14 10 +!ilii 
<13.70 



2 9 +il 



3 5 +1 

—1.0 



1 3 +1 ' 7 
1 - 3 -o.4 



2.1; 



+0. 



1 9 +u 



2 5 +i.o 
^• J -i.o 



2 +1 

^•"-1.0 



[KJ 



183-14J 



262. 
510- 
274 
31 



-161 

-340 
+464 

+2^3 



1+28. 
' 1 -35 

<584 
<610 

244_ 170 

4Q5+398 

324^ 
-tt 

477+54 

-S3 
433 + " 

<607 

<648 

j 34 +229 

223+™ 
324+fe 

<509 
<525 

134_ 63 

148-47 
70 1 +546 

"'-a 
200+g 
<264 
<361 
<446 

209-„, 
371_ 186 

774+341 

-83 

218!l 
<264 
<361 
<446 

134+29 

7Q7+74 
z ''-103 

379™ 

435+ 83 

<408 

7Q1+1075 
z ' 71 -169 

637 +855 
576 +23 ' 
493«§ 
<469 



Velocity relative to z-^s = 2.42630 

2 Best fit values with errors representing the allowed range on 
logW(H2), not the rms error on the fit. Due to the large uncertainties, 
only a range is given for component 1 (see text). 

3 Temperatures corresponding to the best-fit values. In some cases, the 
upper-limit on the ratio ,/V(H 2 , J)/W(H 2 , J = 0) is higher than the max- 
imum expected ratio in case of local thermal equilibrium (see Eq. [T), 
and the corresponding upper-limit on T -j doesn't exist. 



(#7) in these profiles does not correspond to any H2 compo- 
nent and is associated with a weak Ci component. Redshifts 



have been relaxed to allow for the best fit to be obtained. Shifts 
between singly ionised metals components and C 1 components 
are well within errors, however. Voigt-profile parameters given 
in Table H] are well-constrained for species with transitions lo- 
cated outside the Ly-a forest, i.e. Fen, Nin, and Sin. The Sn 
and P 11 transitions are partly blended with intervening Ly-a H 1 
absorptions. Nevertheless, all S 11 components can be constrained 
thanks to the observations of two different transitions. The three 
redder components of Pn are blended, and we cannot derive 
Af(Pii) for these components (see TableHJl. 




-200 



-100 



I''''' 1 " 0.5 



1. ,1 1 ) 1 1 1 1 1 ,1. 1 .1. .1. ) 1 1. 




■ Nl II M74I 
I I I I 1 I I I I I I I I I 1 I I 



SI II M808 
I l_J [__!__! I I ! L_L_ 



' '' '■' 1 I L 



100 -200 -100 

Relative velocity (km s _1 ) 



100 



Fig. 3. Low-ionisation metal transitions in the z a bs = 2.43 DLA 
toward Q 2348-011. Note that the strongest component in the 
metal profiles (at v = km s~') does not show associated H2 ab- 
sorption (see Fig.|4]i. We give upper-limits on A^(P 11) for the three 
redder components of P 11 as they are blended with intervening 
H 1 absorptions. S 11 transitions are also located in the Ly-a for- 
est; but with the help of the two transitions shown above, we can 
still derive column densities or upper-limits for each component. 
Corresponding Voigt-profile parameters are given in Table [4] 



Figure |4] shows the overall absorption profile of highly 
ionised species: C iv and Si iv. It can be seen that the C iv coun- 
terpart of H2 components are generally weak compared to the 
other components. The strongest Civ feature at v = kms -1 
does not correspond to any H2 component but corresponds to 
the strongest component in Fe n. 

The iron-to-sulfur abundance ratio is plotted in Fig. [5] as a 
function of the silicon-to-sulfur abundance ratio. The abundance 
ratios are defined as [X/S]= log^(X)/A?(S) obs -log^(X)/A?(S) . 
A strong correlation between [Fe/S] and [Si/S] is observed, sim- 
ilar to wh at has been observed previously by iPetitiean et al.l 
(2002) and lRodriguez et all ((2006). This correlation is easily in- 
terpreted as a depletion pattern of Fe and Si onto dust grains. 
Note that H2 is detected in components were iron is not strongly 
depleted with [Fe/S] as high as -0.1 in component # 10 and that 
the highest depletions are found in components where H2 is not 
detected, components #5 and 8. However, components #7 and 
# 8 are heavily blended and the distribution of A^(S 11) between 
these two components is uncertain. If we consider these two 
components as a single feature, then the corresponding deple- 
tion would be lower than for component #8. In component #5, 
the S 11 column-density might be overestimated due to poor SNR 
and a small blend (see Fig. [3]). In addition, it can be seen in Fig. 3 
that component # 10 is blended in the wing of component #11. 
Therefore the decomposition is difficult, and the A^(Fe 11) column 
density in this component could be overestimated. 
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Fig. 4. Absorption profiles of Si iv and C iv in the z R bs = 2.43 
DLA toward Q 2348-011. The vertical lines mark the positions 
of low-ionisation metal components in the C iv Si iv and Fe n 
panels. They mark the positions of the seven H2 components in 
the H 2 L0-R0 panel. 



The typical abundance patterns in the Galactic halo, warm 
disk clouds, and cold disk clouds are also indicated in Fig. [5] As 
for Magellanic clouds, the halo pattern represents the depletion 
in individual DLA components best. 



3.4. Chlorine 

When H2 is optically thick, singly-ionised chlorine (CI 11) will 
go through rapi d charge exc hange reaction with H2 to produce 
Hi and Cli(see |juralll974a|). Neutral chlori ne should therefore 
be a good tracer of H? (lJura & Yorkiri978l) and provide addi- 
tional information for constraining the physical parameters of 
the gas. Unfortunately, the wavelength at which the strongest 
chlorine line (CI 1/11347) is expected to fall is not covered by 
our spectrum. It is possible to give an upper-limit on N{Cli), 
thanks to Cli/11088 and Cli/UlOl transitio ns whose oscillator 
streng ths have recently been measured by Son nentrucker et all 
d2006l) . We derive TY(Cli) < 10 13 cirT 2 for each component at 
the 3 cr detection limit. This leads to about 7Y(Cli)//Y(H 2 ) < 10" 5 
for the strongest H2 components (1, 9, and 10). Typical ra- 
tios measured along lines of sight toward the Magellanic clouds 
and the G alactic disk are, re spectively, /V(Cli)/7V(H 2 ) ~ 0.4 - 
3 x 10~ 6 dAndreetal.ll2004l) and 7Y(Cli)//Y(H 2 ) ~ 0.4 x 10" 6 
dSonnentrucker et al.ll2002l) . The non-detection of CI 1 with such 
an upper-limit is therefore not surprising. Singly-ionised chlo- 
rine (CI 11), in turn, has only two transitions at A resl ^ 1063 A and 
1071 A. The corresponding features are blended with H 2 or Ly-or 
lines, but an upper-limit at 3 cr of about N(Cl n) < 10 14 cirT 2 can 
still be derived for components with a relative velocity between 
v -120 kms~' and -30 kms _I (i.e. components 3, 4, and 5). 



Table 4. Results of the Voigt-profile fitting of the singly-ionised 
metal lines. 



# 


Zabs 


Ion (X) 


logN(X) 


b 


Av'x/ci 




v '[kms- 1 ] 






[kms- 1 ] 


[kms- 1 ] 


1 


2.42450(5) 


P11 


12.84± 0.03 


6.6±0.1 


+0.0 




-157 


S11 


14.1 1± 0.02 










Fe 11 


13.73±0.01 










Nin 


12.86± 0.05 










Si 11 


14.30± 0.03 






2 


2.42465(7) 


P11 


12.22± 0.09 


8.1±0.2 


+0.3 




-144 


S11 


14.21±0.02 










Fe 11 


13.70± 0.01 










Nin 


12.92± 0.04 










Si 11 


14.18± 0.04 






3 


2.42494(4) 


P11 


1 1.97± 0.18 


10.7±0.7 


+0.0 




-119 


S11 


< 13.60 










Fe 11 


13.39± 0.01 










Nin 


< 12.50 










Si 11 


13.87± 0.08 






4 


2.42514(4) 


Pn 


12.15±0.09 


5.1 ±0.1 


— 




-101 


S11 


< 13.60 










Fe 11 


13.1 1± 0.01 










Nin 


< 12.50 










Si 11 


13.82± 0.07 






5 


2.42536(8) 


Pn 


11.68± 0.48 


13.5±0.2 


— 




-82 


Sn 


< 14.04 










Fe n 


13.35±0.01 










Nin 


12.81±0.07 










Si 11 


13.85± 0.09 






6 


2.42605(6) 


Pn 


12.32± 0.09 


9.8±0.2 


+0.3 




-21 


Sn 


13.85± 0.05 










Fe 11 


13.34± 0.01 










Nin 


12.71±0.08 










Si 11 


13.93± 0.07 






7 


2.42628(8) 


Pn 


13.09± 0.04 


11.9±0.3 


-0.1 




-1 


Sn 


14.08± 0.07 










Fe n 


14.10± 0.01 










Nin 


13.18± 0.06 










Si 11 


14.55± 0.04 






8 


2.42632(7) 


Pn 


12.45±0.12 


5.7±0.1 


— 




+2 


Sn 


14.49± 0.03 










Fe 11 


13.58± 0.01 










Nin 


<12.50 










Si 11 


14.14± 0.08 






9 


2.42660(0) 


Pn 


12.35± 0.06 


7.4±1.0 


+0.0 




+26 


Sn 


13.78± 0.06 










Fe n 


13.68± 0.01 










Nin 


12.55± 0.09 










Si 11 


14.12± 0.04 






10 


2.42688(8) 


Pn 


12.00± 0.03 


12.0±2.2 


-0.1 




+51 


Sn 


13.90± 0.05 










Fe 11 


14.09± 0.05 










Nin 


12.81±0.07 










Si 11 


14.20± 0.15 






11 


2.42705(4) 


Pn 


12.90± 0.05 


10.5±3.5 


— 




+66 


Sn 


14.04± 0.05 










Fe n 


13.92± 0.04 










Nin 


12.97± 0.05 










Si 11 


14.23± 0.05 






12 


2.42718(3) 


Pn 


blend 


2.8±1.5 


+2.8 




+77 


Sn 


<13.60 










Fe 11 


13.38± 0.01 










Nin 


12.53± 0.08 










Si 11 


13.97± 0.04 






13 


2.42727(9) 


Pn 


blend 


3.4±1.3 


+0.5 




+86 


Sn 


<13.60 










Fe n 


13.47± 0.01 










Nin 


< 12.50 










Si 11 


13.68± 0.08 






14 


2.42741(5) 


Pn 


blend 


6.4±1.2 






+98 


Sn 


<13.60 










Fe n 


13.14± 0.01 










Nin 


< 12.50 










Si 11 


13.86 ±0.06 







1 Relative velocity (in km s 1 ) to the center of the Lyman-o- absorption 
at Zabs = 2.42630. 

Upper-limits are given at 3 cr, except in the case of S n in component 5, 
for which the column-density derived has been considered as an upper- 
limit to take account of a possible blend. 
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Fig. 5. The iron-to-sulfur abundance ratio, [Fe/S], is plotted as a 
function of the silicon-to-sulfur abundance ratio, [Si/S]. The ob- 
served correlation is interpreted as the result of depletion of Fe 
and Si onto dust grains. Different velocity components are num- 
bered from 1 to 14. Filled squares are for components where 
both H2 and C 1 are detected, while filled circles are for compo- 
nents where Ci is detected but not H2. The open circle stands 
for components 7 and 8 considered as a single velocity compo- 
nent. Stars mark the typical depletions observed in the Galactic 
halo clou ds (H) and in cloud s of the Galactic warm disk (W; see 
Table6ofWeltv et al.l lT999). The point corresponding to the de- 
pletion in clouds of the Galactic cold disk would be located off 
of the graph with [Si/S] = -1.3 and [Fe/S] = -2.2. 



4. Determination of physical parameters 

4.1. Excitation temperatures 

The excitation temperature Txy between two rotational levels X 
and Y is defined as 

= H e -Exr/kT X y (1) 
N X gx 

where g x is the statistical weight of rotational level X: g x = 
(23 x + 1 )(2Ix + 1 ) with nuclear spin 1 = for even J (para-Hy and 
1=1 for odd J (ortho-H2), k is the Boltzman constant, and Exy is 
the energy difference between the levels J = X and J = Y. If the 
excitation processes are dominated by collisions, then the pop- 
ulations of the rotational levels follow a Boltzman distribution 
described by a unique excitation temperature for all rotational 
levels. In Fig. [6] we plot, for each H2 component, the natural 
logarithm of the ratios N(U 2 , J = X)/N(U 2 , J = 0) weighted 
by the degeneracy factor gx/go against the energy between ro- 
tational levels. The slope of a straight line starting at the origin 
gives the inverse of the excitation temperatures To-x directly. 

Figure [6] and Table [2] show that 7V 1 is generally much 
lower than the other excitation temperatures. Generally, 7p 1 
is a good indicator of the kinetic temperature dRov et alj r2006; 
iLe Petit et 31.1 12006). However, two effects could act, in opposite 
directions, against this simple assumption. If the J = 1 level is 
not thermalised, then Tqi is a lower-limit on 7\; n . On the con- 
trary, the selective self-shielding of molecular hydrogen in the 
J = 1 level for column densities below NQI2) ~ 10 16 cm" 2 
should make 7oi an upper-limit on 7\i n . Given the large den- 
sities we derive in those components (see §[43}, collisions are 
not negligible and contribute towards thermalising the two lev- 
els. Despite large allowed ranges for individual components, and 




1000 2000 3000 4000 5000 
Ec-x IK] 



Fig. 6. Excitation diagrams for the seven H2 components. The 
grey area shows the range of excitation temperatures that can 
explain the populations of the J > 2 rotational levels. The slope 
of straight lines give directly 1 /To-x- The solid line stands for 
7Vi and dashed lines for the upper and lower limits on 7Vi . 



a relatively large dispersion among the individual components, 
the ratios Af(H 2 , J = 1)/A^(H 2 , J = 0) for the seven H 2 compo- 
nents are all consistent with Tq\ « 140 K, which corresponds 
to the intersection of all the individual Tq\ ranges. This value 
can be taken as representative of the kinetic temperature for the 
derivation of the densities (see § |4.3l l. The kinetic temperatures, 
roughly estimated by Tq\ (see Table d, are nevertheless sim ilar 
to what was found in previous studies lSrianand et al.l d2005l) . us- 
ing the ortho-to-para ratio for thirteen H2 components in seven 
DLAs, found T * 153 + 78 K (see also iLedoux et~alll2006bl 
T ~ 90 - 180 K). These v alues are higher than those mea- 
sured in the ISM (77±17 K; Rachfor detaTI 120021) an d in the 
Magellanic Clouds (82+21 K; iTumlinson et alj|2002h . On the 
other hand, they are similar to what is observed through hi gh- 
latitude Galactic sight lines (124+8 K; iGillmon et"aT1 (120061) or 
ranging from 8 1 K at log N(U 2 ) = 20 to 219 K at log N(U 2 ) = 14; 
Wakker (2006)). This could be due to a high UV radiation field 
but also to different grain properties. Indeed, the photo-electric 
heating of the medium should be more efficient in the case of 
smaller grains. The higher excitation temperatures measured for 
higher rotational levels (J > 2) show that other processes than 
collisions are in play. Rotational levels J = 4 and J = 5 can 
be populated by cascades following UV or formation pumping. 
We indeed assume in the following that these are the main pro- 
cesses populating high-J levels. It is nevertheless also possible 
that the H2 rotational excitation is due to C-shocks or turbulence 
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dJoulainet al.lll998t ICecchi-Pestellini et al. 1 120051) as in diffuse 
clouds. 

4.2. Photo-absorption rate 

We can estimate the total absorption rate of Lyman- and Werner- 
band UV photons by H? in the J = and J = 1 levels, /3q and f}\, 
respectively, following lJural (1 1 9751) . The equilibrium of the J = 4 
and 5 level populations is described by 

P4.0 Po «(H 2> J = 0) + 0. 19 R n(H) n = A 4 ^ 2 n(H 2 , J = 4) (2) 
and 



p 5 ,\ Pi «(H 2> J = 1) + 0.44 R n(H) n = A 5 ^ 3 n(H 2 , J = 5) 



(3) 



where /? is the H2 molecule formation rate, n the proton den- 
sity, p the pumping coefficients and A the spontaneous transition 
probabilities. 

The factor 0.19 (resp. 0.44) means that about 19% ( resp. 
44%) o f H 2 molecules form in the J = 4 (resp. J — 5) level (Jura 
Il974bl) . We use the spontaneous transition probabilities A^ 2 — 
2.8 x 10" 9 s- 1 and A 5 ^ 3 = 9.9 x 10~ 9 s 1 dSpitzeJI 19781) and th e 
pumping efficiencies /J40 = 0.26 and /751 = 0.12 (lJuralll975l) . 
The pumping due to formation of molecules onto dust grains is 
often neglected if densities are lower than 5 x 10 3 cm" 3 , which is 
gener ally the case in damped Lyman-a' systems (Srianand et al. 
2005). However, in the following, we rewrite the above equa- 
tions in a form that only depends on the relative populations of 
H 2 rotational levels. 

The equilibrium between formation and destruction of H 2 
can be written 



fl««(H) = fl diss «(H 2 ) 



(4) 
photo- 



where ^di ss = 0.11 Po, assuming th at 11% of the 
absorptions lead to photo-dissociations d Jurall 1 974bb . 

We c an substitute R n w(H) in Eqs.|2]and[3]and finally obtain 
(see also lCui et aill2005l) : 



Pl.oPo 



jV(H 2 ,J = 0) 
N(H 2 ) 



+ 0.021^o =A4_ 2 



jV(H 2 ,J = 4) 
/V(H 2 ) 



and 



P5,lP 



JV(H 2 ,J= 1) 
N(H 2 ) 



+ 0.049 Pi = A 



5^3- 



jV (H 2 ,J = 5) 
N(M 2 ) 



(5) 



(6) 



One advantage of these expressions is that /?o,i can be determined 
in each component without any assumption to estimate N(H 1) in 
the component. Values obtained for/?o,i in the different compo- 
nents are given in Table 4, where it can be seen that the usual 
assumption Po = Pi is verified. 

We find the photo-absorption rate of H 2 in the Lyman and 
Werner bands p e [7.6 x 10- n ,2.5 x 10~ 9 ] s" 1 (J3 1.1 x 
10" 10 s _1 for the best-fit value) in component 1 ; Pq - lxl0~ 9 s~' 
in component 2; Po =* 1.2 x 10~ 9 s _1 in component 12 and Po ^ 
1.5 x 10 _9 s _1 in component 13. Much lower values are found 
for components with higher H 2 column densities: for component 
9 (Po 1.4 x 10" 11 s^ 1 ), and for component 10 (J3 - 2.4 x 
10 _1 1 s _I ). We also derive an upper-limit on the photo-absorption 
rate for components 3 (/So ^ 2.4 x 10~ 9 s~'). These values are 
spread over about 2 orders of magnitude. 

We plot the photo-abso rption rate, Pq, versus 
logN(H 2 )/iV(Cp in Fig. ISrianand & Petitieanl d 19981) 
suggest that this latter ratio is a good indicator of the physical 
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Fig. 7. Logarithm of the photo-absorption rate as a function of 
log Af(H 2 )/Af(C 1). Squares stand for the best-fit values. Despite 
large allowed ranges on Po and A^(H 2 )/A^(C 1) (dominated by the 
Af(H 2 ) uncertainties), there is a clear anti-correlation between the 
two quantities. 



conditions in the absorber as it increases when the ionisation 
parameter decreases. There is an anti-correlation (R = -0.97) 
between log/?o and log Af(H 2 )/iV(C 1) that is mostly due to 
shielding effects. 

We can correct for the shielding effect (self-shielding and 
dust-shielding) and estimate the UV background intensity out- 
side the cloud. The photo-dissociation rate is related to the am- 
bient UV flux: 



fl diss = 0.11 A, = (4ff) 1.1 x 10 8 7 LW 5 shie id(s" 1 ) 



(7) 



where 7lw (erg s 1 cm 2 Hz 1 sr l ) is the UV intensity at hv 



12.87 eV, averaged ove r the solid angle dAbel et al.l 119971 : 
lHirashita~& Ferrara 2005] ). The term Shield is the c orrection fac- 
tor for H 2 self-shielding dDraine & Bertoldil[T996l) and dust ex- 
tinction. For log NQI2) > 14, it can be expressed as 



1 shield — 



N(H 2 ) 



-0.75 



-o-dMi 



(8) 



The term due to dust extinction, e <TiNi 7 is expressed using 
the optical depth of the dust in the UV, Tuv = 0~dMt: 



t uv = 0.879 



0.1/L/m/ \2gcm 



D 
T(F 



-)( ** ) 

: /\10 21 cm- 2 / 



(9) 



where a is the radius of a grai n, D the dust-to-gas mass ratio, and 
5 the grain material density (Hirashita & Ferrara 2005). Given 
the low H 1 column density, the exponential term is always very 
close to 1, and S shie id - (JV(H 2 )/T0 14 cnr 2 )-°- 75 ; therefore, 



LW 



8 x 10-' l p 



(10) 



* shield 



It can be seen in Table[5]fhat for each component, the allowed 
ranges for 7lw are all much smaller compared to the larger al- 
lowed ranges forySo- This is especially true for components 1 and 
2 (see also Fig. |7]l and caused by errors in N being correlated 
with errors in Doppler parameters. Since absorptions from the 
low J levels are more saturated than those from the higher J lev- 
els, the derived N are more sensitive to a change in b; when b is 
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Table 5. UV radiation field and star-formation rate around indi- 
vidual components 



# 


logjSo 


logjS, 


log y LW 


SsFR 




[s- 1 ] 

n n/:+l.i)6 


[s- 1 ] 


[erg/ s/cm 2 /Hz/sr] 

1 m+0.26 


[M G /yr/kpc 2 ] 

cl 1 n 1 



2 
3 
9 
10 
12 
13 



-9 00 + ^ : ^ 
< -8.62 
-10.85!"° 
-10-62ti? 

_ 8 92+ 0.03 

-8.82 



-n in 



-9.03 + ?f 

< -8.23" 

< -10.67 

< -10.96 

< -9.03 

< -8.74 



' -0.13 

-17 41+ 05 
< -17.74 

-17.82^| 
-17.70^| 



21x10"- 
< 9x10- 
5xl0- 2 
8xl0- 2 
10xl0- : 
13xl0- : 



decreased, log Af(high)/A?(low) decreases, and the derived value 
of /3o gets lower but the shielding correction is larger. Therefore 
the 7 LW values derived for extreme values of b (see § 13.1.21 and 
Table 12) are similar. 

The derived UV intensity is about one to two orders of mag- 
nitude higher than the typical Galactic UV intensity that is about 
7lw.q — 3.2 X 10 -20 ergs _I c m" 2 Hz"" 1 sr"" 1 in t h e sola r vicinity 
dHabinglll968"l) . Using FUSE, iTumlinson et"ai1 d2002b also de- 
rive in the Magellanic clouds UV fluxes 10 to 100 times higher 
than in the Milky Way. This can be interpreted as in-situ star- 
formation close to the absorbing clouds. Note that, after shield- 
ing correction, the scatter in /lw is much smaller than the two 
orders of magnitude variation in /?o from one component to the 
other. It must be recalled that this is derived under the assump- 
tion that high-J excitation is mainly due to UV pumping and 
formation on dust grains. 

The observed ratio of the ambient UV flux to that in the so- 
lar vicinity,^ = Jlw/Jlw,q, is in the range x ~ 30 - 300. This 
should be considered as upper limits, as no other excitation pro- 
cesses than UV pumping and formation on dust is considered 
here. In any case, the observed excitations are still higher than in 
the Galactic halo or the Magellanic clouds. We therefore safely 
estimate that the present UV flux is about one order of magni- 
tude higher than in the ISM. 

Using the relation between the star-formation rate (SFR) 
per unit surface and y, S sfr = 1.7 x \0~^x M yr _1 kpc 2 
(Hirashita & Ferraral 12005), we derive a SFR in the different 
components in the range, Ssfr ~ 5 - 50 x 10 -2 M yr _I kpc~ 2 . 
The detailed values are given in Table It can be seen that 
there is a gradient in the SFR through the DLA profile. The 
derived SFR is higher close to component 1. Note that if the 
star-formation region is larger than about 1 kpc, the associated 
Lyman-q emission sho uld be detectable (e.g M 0ller et alj|2004l 
iHeinmiiller et al.ll2006l) . 



4.3. Ci fine structure 

Neutral carbon is usually seen in components where H2 is de- 
tected. As the ionisation potential of Ci (11.2 eV) is similar to 
the energy of photons in the Werner and Lyman bands, C 1 is usu- 
ally a good tracer of the physical conditions in the molecular gas. 
In the present system, C 1 is detected in all components where H2 
is detected but is also detected in two additional components (# 
6 and 7). It is possible to derive constraints on the particle den- 
sity of the gas from relative populations of the fine-structure lev- 
els of the C 1 ground term (2s 2 2p 2 3 Po,i,2)- In Fig. [8] we compare 
the ratios N(Cf)/N(Ci) and N( C f*j/N(Ci) with mod els com- 
puted with the code POPRATIO dSilva & Viegasll200lh . Taking 
into account the cosmic microwave background (CMB) radia- 



tion with temperature Tcmbr = 2.73 x (1 + z) and assuming 
an electron density of n e = 10~ 4 x «h, kinetic temperatures in 
the range T = 80 - 200 K, and radiative excitation similar to 
that in the Milky Way, we derive hydrogen densities of the or- 
der of n(Hi) a; 100 - 200 cm 3 for the components where H2 
is detected. It is in turn about n(H 1) 50 cirr 3 for components 
without H2 detection if we assume the same kinetic tempera- 
ture. If the temperature is higher in these two components (# 
6 and 7), then densities would be even lower. The results are 
hardly dependent on the exact value of the ambient flux and on 
the electron and H2 densities for this range of particle densities. 
The most influential processes in that case are the radiative ex- 
citation from the CMB and the collisional excitation depending 
on «(H 1) and T. To illustrate this, we also assumed an ambient 
flux 30 times higher than the Galactic flux. The above densities 
for components 6 and 7 are decreased by a factor of about 2 and 
even less for the other components. Note however, that a very 
high UV flux (~ 100 times the Galactic flux) would make the ob- 
served C 1 fine-structure column-density ratios inconsistent with 
the models for the lowest-density components. This again shows 
that the UV flux derived in the previous section is probably an 
upper-limit. 

Assuming temperatures between T ~ 80 K and T ~ 200 K 
also gives very similar results. Using a Galactic ISM tempera- 
ture of T = 80 K leads to a higher inconsistency between the 
N(C f)/N(C 1) and JV(C f*)/N(C 1) ratios. This supports the idea 
already mentioned in the previous section that the kinetic tem- 
perature is higher in this system than in the Galactic ISM (see 
Figs. [6] |8] and Table [6). Note that, because the column densi- 
ties iV(C 1**) are close to the 3 <x detection-limit (log A^(C 1**) = 
12.50) and rather dependent on the continuum placement, den- 
sities derived from iV(C f)/N(C 1) are the most reliable. 

Table 6. Densities derived using models in Fig. [8] with T = 
140 K, and the ratios N(Cf)/N(Ci) (col. 2), AT(C i**)//V(C 1) 
(col. 3). 



# 


JV(Ci*)/AT(Ci) 


n(Hi) [cm- J ] 
N(Cf*)/N(Ci) 


adopted 1 


1 


87-123 


132-162 


105+18 


2 


79-132 


89-151 


106±27 


3 


115-347 


324-776 


231±116 


6 


32-85 


<257 


<59 2 


7 


39-68 


<145 


<54 2 


9 


135-347 


<204 


170±35 


10 


87-123 


79-186 


105+18 


12 


83-138 


162-214 


111+28 


13 


141-240 


309-575 


191+50 



1 Because N(C f) is more reliable than N(C f *), we take the correspond- 
ing n value as the actual n value (col. 4), unless there is an upper limit 
on N(C 1**) that is then taken into account (component 9). 

2 The values for components 6 and 7 are considered as upper-limits, 
because of the higher dependence on the ambient UV flux. 



5. Discussion and conclusions 

We have detected strong molecular hydrogen absorption lines 
associated with the DLA at z abs = 2.4263 toward Q 2348-01 1 in 
seven velocity components spread over ~300 km s -1 and from 
rotational levels J = up to J = 5. We used these observations to 
constrain the physical conditions in the gas. 
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log n(HI) 

Fig. 8. Analysis of the C i fine-structure relative populations us - 
ing results from the POPRATIO code dSilva & Viegasll200ll) . 
Solid and dashed lines give, respectively, N(Cf)/N{Ci) and 
N{Cf*)/N(Ci) versus gas density for a Galactic ambient flux. 
The corresponding thick lines stand for the model with T 
I - 1 K, whereas the thin lines stand for T = 200 K (upper 
curves) and T — 80 K (lower curves). The dashed-dotted lines 
represent the model with an ambient flux 30 times higher than 
the Galactic flux. Boxes show the observed 1 cr range vertically 
for N(C f)/N(C i) and N(C f*)/N(C i) and horizontally the cor- 
responding n values for the model with Galactic ambient flux 
and T - 140 K. They are therefore arbitrarily centred on the 
corresponding curve. Note that the results from different mod- 
els are within errors. The influence of the UV flux is higher for 
components 6 and 7, for which no H2 has been found. The cor- 
responding derived densities are therefore considered as upper 
limits. 



From the fine-structure levels of neutral carbon, we derived 
hydrogen volumic densities of the order of 100 to 200 cm -3 in 
the components where H2 is detected. Densities below 50 cm -3 
are derived in the two components where C 1 is detected b ut not 
H2. This agrees with the results from Srianand et al. (2005). The 
relative populations of the first two rotational levels are all con- 
sistent with a kinetic temperature of about 120 to 160 K, under 
the assumption that Tq\ is a good indicator of 7\j n . The tempera- 
tures derived are higher than typical temperatures in the ISM of 
the Galaxy. Detailed models of diffuse molecular clouds should 
be constructed to understand the origin of this difference. Actual 
kinetic temperatures could even be slightly higher if the J = 1 
level is not thermalised (see, e.g.. lLe Petit et al.l l2006). and then 
densities would also be slightly lower. 

The physical properties of the H2-bearing gas are charac- 
teristic of a cold neutral medium. Molecular hydrogen in this 
DLA traces cold gas but with a higher pressure (of the order of 
p ~ 10 4 ctrT 3 K) than that measured along lines of sight in the 
Galactic ISM. Temperatures and densities are similar to thos e 
observed in the Galactic halo dWakkerl2006tlRichter et al.l2003l) . 

If we assume that H 1 is distributed over the nine C i-bearing 
components proportionally to N(Ci), then the mean Hi col- 



From the relative populations of H2 rotational levels, we de- 
rived photo-absorption rates. Correcting for self-shielding, we 
are able to estimate the incident UV flux on the cloud, which we 
find to be about one order of magnitude higher than in the so- 
lar vicini ty, similar to what is observed through the Magellanic 
clouds (ITumlinson et all2002l) . This shows that there is ongoing 
star formation in this system. 

Searches for the galaxy responsible for the z a b s = 2.4263 
DLA toward Q 2348- 011 have been unsuccessful until now. 
iMannucci et aU (fl998l using the NICMOS3 256 2 MAGIC cam- 
era at the Calar Alto 3.5 m telescope, reported the detection of 
an Ha emission-line galaxy at z - 2.43. The object is located at 
~ 11 arcsec, corresponding to ~ 90 kpc for the standard Q.CDM 
cosmology, from the line of sight toward Q 2348-011 and the 
inferred SFR is 78 M yr _1 . The emission from the correspond- 
ing star-formation region is not larg e enough to explain t he UV 
flux in which the gas is embedded. Bunker et al. (1999) spec- 
troscopically observed Q 2348-011 over the wavelength range 
where possible H a emission at z = 2.43 would be redshifted. 
From the non-detection of any emission, they derived an upper 
limit of about S FR < TAhr 1 M yr _1 for the SFR of any possible 
object exactly aligned with the quasar. However, their observa- 
tions are not very sensiti ve (~10 -15 erg/s/cm 2 ). More recently 
IColbert & MaTkanl (120021) obtained coronographic H-band im- 
ages of the field with the Near-Infrared Camera and Multi-Object 
Spectrograph mounted on the Hubble Space Telescope. With a 
5cr detection limit H = 22, they detected three objects at 3.3, 9.7, 
and 10.7 arcsec or 19, 55, and 61 h' 1 kpc from Q 2348-011. If 
at the same redshift as the DLA, these objects are about 1 to 2 
magnitudes brighter than an L, galaxy. In particular the authors 
reject the object closer to the line of sight as a possible DLA 
galaxy candidate because of its brightness. 

The origin of the UV excess in the vicinity of the gas is there- 
fore unclear. Either one or several objects within 1 1 arcsec from 
the quasar are indeed at the redshift of the DLA and provide the 
inferred UV flux through star-formation activity or the region of 
star-formation is located on the line of sight to the quasar that is 
exactly aligned with it. In that case, Lyman-a emission should 
be observed at the bottom of the DLA trough. The SDSS spec- 
trum of Q 2348-01 1 is too noisy to claim anything in this matter. 
It will thus be very interesting (i) to determine the spectroscopic 
redshifts of the objects that are located close to the line of sight 
and (ii) to perform deep intermediate-resolution spectroscopy of 
the quasar to search for any Lyman-o- emission. 
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umn density is about 10 cm in individual components. 
The molecular clouds must therefore have an extension along 
the line of sight that is smaller than 1 pc. This supports the 
view that H2-bearing gas is clumped into small clouds (see, e.g., 
iHirashita et alll2003b . 
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Fig. 9. H2 rotational level J = 0. The centre of the velocity scale 
is set at z a b s = 2.4263. Continuum is seen depressed for H 2 L2- 
0R0 and H2 L1-0R0 because of the z a bs = 2.615 Ly-/3 absorp- 
tion. Components are marked by dashed lines. The strong ab- 
sorptions at v ~ -10 - km s _1 is due to J = 1 transition lines. 
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Fig. 12. H 2 rotational level J = 3. H 2 L1-0R3 is blended with 
FeiLil096. 
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Relative velocity (km s *) 

Fig. 10. H2 rotational level J = 1 . Continuum is depressed in 
panels H 2 L2-0R1 to H 2 L1-0R1 because of z abs = 2.615 Ly-/3 
absorption. The J = 1 components are marked by dashed lines, 
the other ones come from other transitions (in general J = - 2). 
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Fig. 13. H 2 rotational level J = 4. Continuum is depressed in 
some panels by z a bs - 2.6 Ly-yS absorption. The feature in the 
H 2 Ll-0 P4 panel at v > -50 km s might be due to Civ absorp- 
tion at z a b s - 1.45. 
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Fig. 11. H 2 rotational level J = 2. H 2 L1-0P2 is blended with 
Fen/11096. 
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Fig. 14. H 2 rotational level J = 5. Most lines are blended with 
Ly-a or with z a b s - 1.45 Civ absorptions. Column densities are 
mainly constrained by H 2 L3-0 R5 for components 1 & 2 and by 
H 2 L4-0P5 (resp. H 2 L2-0R5) for component 9 (resp. 10). 
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ABSTRACT 

Aims. We aim at deriving the physical conditions in the neutral gas associated with damped Lyman-a- systems using observation and 
analysis of IF and C i absorptions. 

Methods. We obtained a high-resolution VLT-UVES spectrum of the quasar Q 2348-01 1 over a wavelength range covering the metal 
and molecular absorptions from the log N(H i) = 20.50+0.10 damped Lyman-a system at z a bs = 2.4263. We detect IF in this system 
and measure column densities of IF, Ci, Ci*, Ci", Pn, Sn, Fen, Sin and Nin. From the column density ratios and in particular 
the relative populations of the IF rotational and C i fine-structure levels, we derive the physical conditions in the gas (abundances, 
dust-depletion, particle density, kinetic temperature and ionizing flux). 

Results. Molecular hydrogen is detected in seven components in the first four rotational levels (J = 0-3). Absorption lines of IF 
J = 4 and 5 rotational levels are detected in four of these components. This leads to a total molecular fraction of log / = — 1.69^5 5L 
Fourteen components are needed to reproduce the metal-line profiles. The overall metallicity is found to be -0.80, -0.60, -1.17, 
-0.89, -0.45+0.10 for, respectively, [Si/H], [S/H], [Fe/H], [Ni/H] and [P/H]. We confirm the findings by others that there is a clear 
correlation between log N(¥e n)/N(S n) and log./V(Siii)/W(S n) from different components indicative of a dust-depletion pattern. The 
dust depletion in the H 2 components is not large and ranges from log iV(Fe n)/N(S n) = -0.8 to -0.2 when the most important depletion 
(up to log 2V(Fe u)/N(S n) = —1.2) occurs in components where IF is not detected. The gas in IF bearing components is found to be 
cold but still hotter than similar gas of our Galaxy (T ~ 140 K instead of typically 70 K), dense (n ~ 50 - 200 cm -2 ), and embedded 
into a moderate to high UV field, about one to two orders of magnitude higher than in the solar vicinity. This suggests the gas in 
IF. -bearing DFAs is clumpy, and star-formation probably occurs in the associated object. 
Conclusions. 

Key words, galaxies: ISM - quasars: absorption lines - quasars: individuals: Q 2348-01 1 



1. Introduction 

High-redshift damped Lyman-a systems (DLAs) detected in ab- 
sorption in QSO spectra are characterized by their large neutral 
hydrogen column densities, N(Hi) > 2 x 10 20 cirT 2 , similar to 
what is measured through local spiral disks (?). The correspond- 
ing absorbing clouds are believed to be the reservoir of neutral 
hydrogen in the Universe (e.g. ?). Though observational studies 
of DLAs have been pursued over more than two decades, im- 
portant questions are still unanswered, such as (i) the amount of 
in-situ star-formation activity in DLAs, (ii) the connection be- 
tween observed abundance ratios and the dust content, and (iii) 
how severe is the bias due to dust obscuration in current DLA 
samples. 

One way to tackle these questions is to derive the physical 
conditions of the gas. This is possible because, thanks to the shift 

* Based on observations carried out at the European Southern 
Observatory (ESO) under prog. ID No.072.A-0346 with the UVES 
spectrograph installed at the Very Farge Telescope (VLT) on Cerro 
Paranal, Chile. 



of the lines towards the red part of the spectrum, most of the 
species can be observed at once in the same spectrum. In the case 
of the Galactic interstellar medium (ISM), very few lines of sight 
are available with the information from molecular hydrogen and 
metal species together. 

Although molecular hydrogen is the most abundant molecule 
in the Universe, it is very difficult to detect directly however. The 
only way until now comes from the use of the UV Lyman and 
Werner absorption signatures left in the spectra of background 
sources such as quasars (?) or, more recently y-ray burst after- 
glows (?). 

Formation of molecular hydrogen is expected on the surface 
of dust grains when the gas is cool, dense and mostly neu- 
tral. At Hi column densities as large as those measured in 
DLAs (logV(Hi) > 19.5), H2 molecules are conspicuous in 
our Galaxy: gaseous clouds with log N{Hi) > 21 usually have 
log ATH2) > 19 (see ??). The amount of H2 molecules seen in 
DLAs is much less than this. ? find that about 15% of DLAs 
show detectable H2 with molecular fraction in the range -4 < 
log / < — 1. The fraction increases with metallicity and reaches 
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50% for [X/H] > -0.7 (?). The fact that H 2 is less conspicuous at 
high redshift is easily explained as a consequence of lower metal 
content and therefore lower dust content together with higher 
ambient UV flux. 

Detailed analysis of H2-bearing DLAs have been pursued 
in some cases. ?, using R ~ 30000 STIS data from the HST, 
recently derived a low hydrogen density « H ~ 0.2 crrT 3 and 
a electron temperature of T e « 140 K in the H 2 absorber at 
Zabs = 1-7765 toward Q 1331 + 170. Other authors however de- 
rive higher densities. By studying the spectrum of the z a bs - 2.81 
DLA toward PKS 0528-250, ? derived a kinetic temperature 
of the cloud T ~ 200 K and a density n ~ 1000 crrT 3 . They 
suggested that the ratio N(H2)/N(Ci) could be a good indica- 
tor of the physical conditions of the absorber. ? derived gas 
pressure, particle density, excitation temperatures and UV am- 
bient fluxes for four H2 components at z a b s - 1-973 toward 
Q 0013-004. They show that whenever H2 is detected, the par- 
ticle density is high («h ~ 30 - 200 cm -3 ), and the temperature 
is low (T ~ 100 K). In addition, the UV ambient radiation field 
is highly inhomogeneous. ?, using 33 systems from the high- 
resolution (R ~ 45 000) UVES survey from ? - among which 8 
show detectable H2 absorption lines - derived the mean kinetic 
temperatures, the density of hydrogen and the electronic density 
together with the UV ambient flux in DLAs. They conclude the 
mean kinetic temperature in these clouds is higher (~ 150 K) 
than in the local ISM (~ 80 K), the UV flux is of the order 
or higher than the mean UV flux in the ISM of our galaxy, 
and that the densities (about 10-200 cirT 3 ) and the ionisation 
rate are similar to the one in the Galactic cold neutral medium. 
More recently, gas temperature, T ~ 90 - 180 K, and density, 
«h < 50 crrT 3 , were derived in a DLA at z abs 4.224 which 
corresponds to the highest redshift at which this kind of study 
has been done up to now (?). 

We recently reported the detection of H2 in the z a b S = 2.4263 
DLA toward Q 2348-011 (?) and here we present the detailed 
analysis. This system has one of the highest molecular frac- 
tions ever measured at high redshift, with several components 
in which physical conditions can be derived. The paper is organ- 
ised as follows. In Section 2, we describe observations and data 
reduction. In Section 3, we give details on the data and present 
the fits to the H2, C 1 and low-ionisation metals absorption lines. 
In Section 4, we discuss the relative populations of the H2 ro- 
tational levels, and the fine-structure levels of Ci, deriving in- 
formation on the physical state of the gas. Finally, we draw our 
conclusions in Section 5. 

2. Observations 

The quasar Q 2348 -011 was observed with the Ultraviolet and 
Visible Echelle Spectrograph (UVES, ?) mounted on the ESO 
Kueyen VLT-UT2 8.2 m telescope on Cerro Paranal in Chile 
in visitor mode on October 29 and 30, th , 2004, under good 
seeing conditions (FWHM < 0.8 arcsec) and at low airmass 
(AM < 1.2). Three 5400 s exposures were taken with Dichroic 
#1 simultaneously with the blue and red arms. We used a 1 .0 arc- 
sec slit width and a pixel binning of 2 x 2, leading to a spectral 
resolution R ^ 47 000. The data were reduced using the UVES 
pipeline (?) which is available as a context of the ESO MIDAS 
data reduction system. The main characteristics of the pipeline 
are to perform a precise inter-order background subtraction, es- 
pecially for master flat-fields, and to allow for an optimal extrac- 
tion of the object signal rejecting cosmic rays and performing 
sky subtraction at the same time. The pipeline products were 
checked step by step. 



The wavelength scale of the reduced spectrum was then con- 
verted to vacuum-heliocentric values and the portions of the 
spectrum corresponding to different settings were each rebinned 
to a constant wavelength step. No further rebin was performed 
during the analysis of the whole spectrum. Individual 1-D ex- 
posures were scaled, weighted and combined together using the 
scombine task of the noao onedspec package of IRAF. Standard 
Voigt-profile fitting methods were used for the analysis to deter- 
mine column densities using the oscillator strengths compiled in 
Table 1 of ? for metal ions and the oscillator strengths given by 
? for H 2 . 

We adopted the Solar abundances given by ? based on the 
meteoritic data from ?. H 2 transitions wavelengths are from ?. 

3. The DLA at z abs - 2.4263 toward Q 2348-011 

There are two DLAs at z a b s = 2.42630 and z a bs = 2.61473 to- 
ward Q 2348-011. These DLAs were discovered as strong Ly- 
a absorptions in a low spectral resolution (~ 5 A) spectrum 
by ?, who also observed strong metal lines associated to the 
z a bs = 2.4263 DLA. From Voigt-profile fitting to the damped 
Lyman-a line, we find that the total neutral column density of 
the system at z abs = 2.4263 is logJV(Hi) = 20.50 + 0.10. A si- 
multaneous Voigt-profile fitting of the Lyman-a and Lyman-/? 
lines has also been performed for the system at z a bs = 2.61473. 
Note that the z a b s = 2.61473 Lyman-/? line falls in the wave- 
length range where H 2 lines at z abs = 2.4263 are detected, see 
Fig. ??. Total neutral column density for this system is log 7V(H i) 
= 21 .30+0.08. No molecular hydrogen is detected at this redshift 
down to log JV(H 2 ) < 13.7, or log/ < -7.2. 

3.1. Molecular content 

H2 lines at z a b s = 2.4263 are strong and numerous (see Fig. ??). 
Although they are all redshifted in the Lyman-a forest, the spec- 
tral resolution of our data (R ^ 47000) is high enough to allow 
unambiguous detection and accurate determination of the line 
parameters. 

3.1.1. Overall Fit 

The continuum normalisation has been performed using first or- 
der broken lines over wavelength ranges of about 50 to 100 A 
between A =* 3600 A and A 3850 A, using a routine based on 
the IRAF task continuum. Note that only the Lyman band is cov- 
ered by our spectrum. 

The absorption profiles of the once-ionised metal species 
have been fitted using fourteen components and we use their 
numbering to refer in the following to any component of the 
cloud. Seven of these components, spread over about 300 km s _1 , 
show associated H2 absorption. Multiple component Voigt- 
profile fittings have been performed simultaneously on all H2 
absorption lines seen redward of the Lyman break, also includ- 
ing the Voigt profile of the Ly-/3 line at z a b s = 2.61473. 

The solution of a multi-component Voigt-profile fit is often 
not unique, especially when absorption lines are blended and 
sometimes saturated. In that case, it is difficult to constrain both 
the Doppler parameters and the column densities at the same 
time for all species. We therefore obtained a first guess of the 
molecular parameters by fitting the lines of the four rotational 
levels, J = to J = 3, using Doppler parameters and redshifts 
derived from the C 1 profiles which are not saturated and located 
outside the Lyman-a forest. We then froze the corresponding pa- 
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Fig.l. Portion of the normalised spectrum of Q2348-011. The spectral resolution is about R = 47000. The H2 spectrum derived from Voigt 
profile fitting is superimposed to the quasar spectrum. The strong feature at A & 3707 A is the Ly-/? absorption line from the DLA at z a bs = 2.6147. 
The two similar features at A ~ 3784A and A ~ 3790A are likely to be Civ absorptions at z abs ^ 1.45. H 2 components are numbered for the J = 
rotational level. Transition names are abbreviated: H 2 Lx-OYz is written LxYz. The fit parameters are given in Table ??. 



rameters and added the J = 4 and J = 5 features. After a satis- 
factory fit to all H2 transitions was found in these conditions, we 
relaxed all parameters one by one, until a good fit was obtained. 
Then, a range of column densities is derived by fitting the ab- 
sorption features with the lowest possible Doppler parameters 
and the highest ones. The errors given in column 4 of Table ?? 
therefore correspond to the ranges of column densities for which 
the fit is acceptable. 

The quality of the overall fit can be estimated on Fig. ?? and 
from Fig. ?? to Fig. ??. Results are presented in Table ??. The 
total H2 column density integrated over all rotational levels is 
log/Y(H 2 ) = 18.52+9??, corresponding to a molecular fraction 



log / = -1.69+HZ. We also derive an upper limit on the column 
density of molecular HD, leading to log /Y(HD)//y(H 2 ) < -3.3. 
We describe below the individual components. 

3.1.2. Comments on individual components 

Components 1 & 2, z a bs = 2.424497 and Zabs = 2.424638: These 
two components are heavily blended and appear as a single 
saturated feature for low H2 rotational levels, but are better 
separated for J > 4, in C 1 and low-ionisation metal lines (See 
Fig. ??, ?? and ??). Doppler parameters, b, are mainly con- 
strained by the H 2 L0-0R0, H 2 L1-0P2, H 2 L1-0P3 profiles, 
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and the transitions from the J > 4 rotational levels, but remains 
uncertain. The first guess of the H2 profile using redshifts from 
Ci, Ci* and Ci** transitions was extremely helpful in this case. 
It must be emphasize that the blends of the two components 
can be fitted using either low b values, leading to high column 
densities or higher b values, leading to lower column densities. 
No damping wings are present which would be in favour of 
very low b values. This is why we give a column density range 
for component 1. The best fit however is obtained with b = 2.9 
kms -1 and column densities log/V(H 2 , J) = 17.15; 17.70; 
17.00; 17.60; 15.40 and 15.10 for J = to J = 5, respectively. 
These are the values we use in the following. Rotational levels 
are unambiguously detected up to J = 5. Possible absorptions 
from J = 6 and J = 7 levels, unfortunately, are always redshifted 
on top of strong H 1 Lyman-a features or in regions of bad SNR. 
Upper limits are given in Table ??. 

Component 3, z a b s = 2.424917: This component is the weakest 
in the system, and is not saturated. Doppler parameter is 
therefore mainly constrained by the strongest transition lines 
that are not blended: H 2 L2-0R0, H 2 L1-0R0 and H 2 L3-0P1. 
Only the first 4 rotational levels are detected. Components 9 & 
10, Zabs = 2.426591 and z a bs = 2.426885 : These components 
have high H 2 column densities. They are well separated and 
the Doppler parameters are well constrained by the width of 
the isolated and non-saturated absorption features: H 2 L0-0RO, 
H 2 L0-0R1 and most of the P-branch (AJ = -1) Lyman-band 
transitions from the J = 2 to J = 4 rotational levels. Note 
that component 9 of a particular transition is often blended 
with components 1 and 2 of another transition from the next 
rotational level (see e.g. H 2 L2-0R0 in Fig. ??). The six first 
rotational levels (J = to J = 5) are detected. Column densities 
in the J = 5 rotational level are not well constrained (see 
transitions H 2 L4-0P5 and H 2 L2-0R5 on Fig. ??) in these 
two components. Same procedure has been applied as for 
components 1 and 2 for J = 6 and J = 7 levels: upper-limits have 
been derived as the largest value consistent with the observed 
spectrum. 

Components 11 & 12, z a bs = 2.427135 and z a bs = 2.427293 : 
Absorptions from rotational levels J = to J = 3 are detected 
in the two components. J = 4 and J = 5 rotational levels are 
detected but the corresponding column densities are uncertain 
and one should probably consider these values as upper limits 
only (See Table ??). 



3.2. Neutral carbon 

Nine components are detected in the complex and sometimes 
highly blended profiles of neutral carbon. Lines from the 
fine-structure triplet levels of the ground state (2s 2 2p 2 3 Po,i,2) 
are detected. In Fig. ?? we show the fit to the Ci absorp- 
tion features around A « 5346 A and A « 5676 A with 
contributions of the eleven transitions: C 1/11560, C 1/11656; 
C r*/U560a, C i*/U560b, C i*A1656, C 1* /11657a, C 1* /11657b and 
Ci"/U561a,Ci"/U561b, Ci**A1657, Cf"A1658. We first per- 
formed a fit with seven Ci components corresponding to the 
seven H 2 detected components. We left however the exact posi- 
tions of the C 1 components free. We then added C 1* and finally 
C 1** transitions. To reproduce the overall profile well we had to 
introduce two weak C 1 components (# 6 and 7) that are not de- 
tected in H 2 , but are detected in the low-ionisation metals. We 
then then re-performed a fit with nine components, in the same 
order: first Ci only, then adding Ci* and finally Ci**. Finally 



Table 1. Voigt-profile fitting results for different rotational levels of the 
vibrational ground-state transition lines of H 2 . 
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1 Velocity relative to Zubs = 2.42630 

2 Best fit values with errors representing the allowed range on 
log A^I-b), not the rms error on the fit. Due to the large uncertainties, 
only a range is given for component 1 (see text). 

3 Temperatures corresponding to the best-fit values. In some cases, the 
upper-limit on the ratio Af(H 2 , J)/7V(H 2 , J = 0) is higher than the maxi- 
mum expected ratio in case of local thermal equilibrium (see Eq. ??), 
and the corresponding upper-limit on T -j doesn't exist. 



all transitions are fitted together. Small shifts (see Tab. ??) in 
velocity space between C 1 and H 2 components are within uncer- 



